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•II - SOLUTION BY GREEN'S FUNCTION 
By John D. Stanitz 


SUMMARY 

Methods of solution by Green's function. are developed for the design 
of two-dimensional unbranched channels with prescribed velocities as a 
function of arc length along the channel walls. The methods apply to 
incompressible and linearized compressible, nonviscous irrotational flow. 
One numerical exaii 5 )le is presented for an accelerating elbow with 
linearized compressible flow. The elbow shape obtained from the solution 
by Green’s function is the same as that .obtained from a solution by 
relaxation methods for the same prescribed conditions ^ The time required 
for the calculations is considerably less for solutions by Green's 
function. 


INTRODUCTION 

In this report a general method of design is developed for two- 
dimensional, compressible or incompressible, nonviscous irrotational 
flow in unbranched channels with prescribed velocities as a function of 
arc length along the channel walls. The design of channels with pre- 
scribed velocities is important because: (i) boundary- layer separation 

losses can be avoided by prescribed velocities that do not decelerate 
rapidly enough to cause separation, (2) shock losses in compressible 
flow and cavitation in incompressible flow can be avoided by prescribed 
velocities that do not exceed certain ■ maximum values dictated by these 
phenomena, and (3), for compressible flow the desired flow rate can be 
assured by prescribed velocities .that do .not result in "choke flow" . 
conditions. 

In Part I of this report (reference 1) solutions were obtained by - 
relaxation methods. This method of solution results in complete infor- 
mation concerning the distribution of flow .conditions throughout the 
channel and can be used to obtain solutions for incompressible flow 
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for two types of compressible flow: the general type with arbitrary- 

value for the ratio of specific heats y (1.4, for exanqjle) and the 
linearized type with p equal to -1.0. 

In the present report solutions are obtained by Green's function. 
This method of solution is limited to incompressible and linearized 
(p = -l.O) compressible flow, but the method is more rapid than relax- 
ation methods, provided information within the channel is not required. 
TSie method of solution is developed fbr the channel walls only although 
the method can be extended to determine the shape of streainlines within 
the channel. 

The design method reported herein -was developed at the NACA Lewis 
laboratory during 1950 and is part of a doctoral thesis conducted with 
the advice of Professor Ascher H. Shapiro of the Massachusetts Institute 
of Technology. 


METHOD OF SOLUTION 

The design method is developed for two-dimensional channels with 
prescribed velocities along the channel walls. The prescribed velocity 
is arbitrary except that stagnation points (zero velocity) cannot be 
prescribed. This exception limits the design method to unbranched 
channels. In the present report the method of solution is by Green's 
function in conjunction with a formula derived (elsewhere) from Green's 
theorem. 


Preliminary Considerations 

Assumptions. - The fluid is assumed to be nonviscous and either 
compressible or inconpressible. ,If the fluid is conpressible, the 
ratio of specific heats p is assumed to be -1.0, so that the differ- 
ential equations describing the flow are linear. The flow is assumed to 
be two-dimensional and irrotational. 

Physical plane. - The flow field of the two-dimensional channel is 
considered to lie in the physical xy-plane where x and y are 
Cartesian coordinates expressed as ratios of a characteristic length 
equal to the constant channel width do-vmstream at infinity. (All symbols 
are defined in appendix A.) 

At each' point in the channel the velocity vector (fig. 1) has a 
magnitude Q and a direction 6 where Q is the fluid velocity 
expressed as the ratio of a characteristic velocity equal to the constant 
channel velocity downstream at infinity. For compressible flow, the 
velocity Q is related to the velocity q by 
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q (1) 

where q is the velocity expressed as a ratio of the stagnation speed 
of sound and the subscript d refers to conditions downstreeim at 
infinity. 

Stream function f . - If the condition of continuity is satisfied, 
a stream function ijr can be defined such that for incompressible flow 

■ t 

d¥ = I di|f (2a) 


where, from Part I, 


,d\lf = Q dn 


(2b) 


where n is distance in the xy-plane measured normal to the streamline 
and expressed as a ratio of the channel . width downstream at infinity. 
For linearized compressible flow (y = -1.0) 


d¥ - - ^ 

* 2 Aip 


where, from Part I, 


(2c) 


dijf* = p*q* dn (2d) 

and where Ailr* is the value of \|f * along the left channel wall when 
faced in the direction of flow if the value of ilr* along the right wall 
is arbitrarily equal to zero. The value of Ai|f* is obtained by inte- 
grating equation (2d) across the channel at a position far downstream 
where flow conditions are uniform 

At* = P^*%* (2e) 

From part I, p* is related to the density p., expressed as the ratio 
of a characteristic density equal to the stagnation density, by 
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where 



in which 



(2h) 


(2i) 


where the subscripts a and b refer to quantities related to any two 
selected values of velocity {q and q, , respectively) for which veloc- 

ities the densities given by equations (2h) and (2i) are equal to the 
densities p given by' equations (2f) and (2g) , Also, from Pajj-t I, 
q* is related to q by 

q* = (2j) 


where 


k2 



(2k) 


For each prescribed velocity distribution along the' channel walls there 
are an infinite number of linearized compressible flow solutions, 
depending on the selected values of qg^ and q^ in equations (2g) and 

(2k) . However, for values of q^ and q-j^ within the range of q pre- 
scribed along the channel walls (and therefore everywhere in the channel), 
the solutions, that is, channel shapes, probably differ only in small 
detail (Part I) . 
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The values of and q^ might, for exaii 5 )le, be selected to equal 

the maximum and minimum values of q (which values of q must occur on 
the channel walls and are therefore known) . Also, the values of qg^ 
and q-jj might be selected to equal the upstream and downstream veloc- 
ities q.^ and q^. In this case the upstream and downstream channel 
widths would then satisfy continuity for a gas- with the correct 
(arbitrary) value of y (1.4, for example) . If the upstream and down- 
stream velocities are equal, their value and the value of some other 
velocity (the maximum or minimum velocity, for example) can be selected 
for qg and q-j^ or, if desired, q^ can be equal to q-j^ so that 


q 


a 


q 


and 


and, from Part I, 


and 


Pa = Pb = P 



2 

q 



( 21 ) 


(2m) 


Equations (2a) and (2c) define the stream function 'W for incom- 
pressible and linearized compressible flow, respectively. For both 
types of flow W varies from zero along the right side of the channel, 

when faced in the direction of flow, to along the left side of the 
channel . 

Velocity potential $ , - If. the condition of irrotational fluid 
motion is satisfied, a velocity potential $ can be defined such that 
for incompressible flow 
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where, from Part I, 


dCp = Q ds (3b) 

where s is distance in the xy-plane measured along the streamlines 
and expressed as the ratio of channel width downstream at infinity. 

For linearized compressible floy 


n dcp” 

'!«’= sSp 


(3c) 


where, from Part I, 


dCp* = q* ds - (3d) 

s 

Equations (3a) and (3c) define the velocity potential <I> for incom- 
pressible and linearized con^ressible flow, respectively. 


Outline of design method . - Solutions for two-dimensional flow are 
boundary- value problems. That is, the solutions depend on known con- 
ditions imposed along the boundaries of the problem. In the inverse 
problem of channel design the geometry of the channel walls in the 
physical xy-plane is unknown. This unknown geometry apparently precludes 
the possibility of solving the problem in the physical plane and neces- 
sitates the use of some new set of coordinates, that is, a transformed 
plane, in which to solve the problem. These new coordinates must be 
such that the geometric boundaries along which the velocities are pre- 
scribed are known in the transformed plane. It is also necessary, for 
the method of solution employed in this report, that the coordinate 
system of the transformed plane be orthogonal in the physical plane. 

A set of coordinates that satisfies these requirements is provided by 
$ and ¥, which are orthogonal in the physical xy-plane and for which 
the geometric boundaries are known constant values of ¥ (equal to 0 

and ^) in the transformed ^>¥-plsne. The distribution of velocity as 

a function of $ along' these boundaries of constant ¥ is known because, 
if 


Q = Q(s) 
or 

q = q(s) 

is prescribed, equation (3a) or (3c) integrates to give 


= $(s) 
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from which 

Q = Q($) 


or 


q = q($) 

The technique of the proposed channel- design method is therefore to 
solve for the physical x,y- coordinates of the channel walls in the 
transformed $Tir-plane where the prescribed boundary conditions for the 
two-dimensional, flow problem are known. 

Channel wall coordinates . - From Part I the distribution of channel 
wall coordinates x and y along the boxmdaries of constant W equal - 

to 0 and ^ in the transformed ^ijf-plane is given by 


and 



(4a) 



for linearized compressible flow, and for incompressible flow 


(4b) 


and 


^ _ 2 I cos 9 
n Q 


d® 


2 sin 9 ^ 

jt q 


(5a) 


(5b) 


where the constants of integration are selected to give known (specified) 
values of x or y at one value of $ along each boundary. Because 
q* and Q are known functions of^$ from the prescribed velocity as a 
function of arc length along the channel walls, the shape of the channel 
walls in the physical xy -plane is given by equation (4) or (5) if 9 is 
determined as a function of $ along the channel walls (f equals 0 and 

. In this report ,the solution for 0 as a function of, $ along the 
channel walls in the -plane is obtained by Green’s function. 
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Solution by Green's Function 


Continuity . - From Part I the continuity equation becomes in the 
transformed $^plane 


+ 3^= 0 




( 6 a) 


where for ihcon5)ressible flow 

V = Q. 

and for linearized compressible flow 


( 6 b) 


V = 


1 + + q*2 


( 6 c) 


Irrotational motion . - From Part' I the equation for irrotational 
motion becomes in the transformed OW-plane 

S logg V ^0 ^ 

5 ? ^ = ° 

Integral eq-uation for 0 ('I’o^^o)* “ equations ( 6 a) and ( 7 ) 




( 8 ) 


so that from appendix B the value of 6 at a point ($q,,¥^) within, 

or on, the channel walls in the transformed $lir-plane is given by the 
integral equation 


e($o,¥o) 


f 

/ a logg y\ 

) « 1 

^ a log 

r . ^ 

^0 






I 

-00 

L\ / 

2 

\ / 

- 


d$ 


( 9 ) 
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where the subscripts 


0 and — 
2 


refer to the channel wall boimoaries 


along' which W is 0 and respectively, and G is the Green's 
function of the second kind for the channel, which is an infinite strip 
of width ^ extending in the $-direction to ±». 


Green's function G . - The Green's function of the second kind G 
for the infinite channel in the $W-plane is given along 'the channel 

wall boundaries (W equals 0 and •^) by (appendix C) 


or rt = -loge 


cosh^ (^-^o) ” cos^ 


(10). 


where ($,¥) is any point on the channel wall boundary and 1^ 

the point in the channel or on the boundary at which 9 is- to be deter- 
mined . 


Numerical integration for - From equations (9) and (lO) 


‘2jt e($ ) 

' o o' 



cosh^. 


sin^ ¥ 
o 



d($-$ ) - 



cosh 




( 11 ) 


in which the independent variable of integration has been changed from 
d^> to d($-<I>Q) so that the origin, for purposes of integration, lies 
at rather than $ = 0. If for small changes in (<i>-$Q), that is 

S logg V ' 

for small A$, the term may be considered constant and equal 

to its average value over the interval A<5, then 


S logg V A logg V 

— ^ — = — aT~ 


and equation (11) becomes 
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00 

z 


($-$o)=-Oot 


[A logg V 

"a¥“ 






log„ cosh^($-$ ) - sin^ ¥ d($-$)> 

L °J ° J n 


^ (2!i log V r 

4-. { A^~ / logg cosh^ ($-(!> J 


- cos*- 


^o] 


( 12 ) 


where the summation sign is understood to mean that the quantity within 
the braces is summed over the entire range of (0-$^) between ±<». 


Equation (12) determines 9 at any point in the flow field 
(channel). For a point (<I>^,'W^) on the channel walls is equal 

to 0 0^ ^ integrands in equation (12) become 


2 logg cosh |($-$o)l 
or 


2 logg sinh |($-$g)| 
so that e’quation (12) becomes 


\ 


OO 

2 : 






/A logg V 
A^^ 



(13a) 


where 

\ 




(13b) 
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where 


= a if Vo =■ 0 

2 

In = P if 

2 

lo = a If = I 

Iq = P if ''o = ° 


KJ-Jg)! 


a = 


logg cosh|($-$o)| d|($-'l'o)( 



logg sinh|($-$^)| d|(<I>-$^)| 


(13c) 


.(13d) 


(I3e) 


where the + signs apply for positive values of ($-^q) and the - signs 
apply for negative values of (<I)-$o) . Methods of evaluating a and 3 
are given in appendix D and tabulated values are given for a wide range 
of K*^“'^*o)l table I. Equation (l3a) determines 6(^o^^) any 
point on the channel wall boundaries. Thus from equations (4a) and (4b) 
or (5a) and (5b) the coordinates for the channel wall shape in the 
physical xy-plane can be determined. 


NUMERICAL PROCEDURE 

The numerical procedure for the channel design solution by Green's 
function is the same, except for minor details, for incompressible and 
linearized compressible flow. The stepwise procedure is outlined as 
follows : 

(1) For incompressible flow the velocity Q and for linearized 
compressible flow the velocity q, or which is the same thing the 
velocity Q and the constant downstream velocity q^, are specified 
as functions of arc length along the channel walls 


Q = Q(s) 


(14a) 
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or 


q = q(s) ' (14b) 

where s is arbitrarily equal to 0 at that point along one channel 
wall where the velocity first begins to vary. , 

(2) Compute V as a function of s from equations (6b) and (I4a) 
for incori 5 )ressible flow or from equations (2 j) , (2k) , (6c) , and (I4b) 
for linearized compressible flow. 

V = V(s) ' (15) 

(3) Compute $ as a function of s from equations (3a) and (3b) 

for incompressible flow or from equations (2e) , (3c) / (3d ) , and (l4b) 
for linearized compressible flow. In equation (2e) p^* is obtained 

from equations (2f) to (2i) . For arbitrary distributions of Q or q 
equation (3a) or (3c) is integrated numerically using, for example, 
Simpson’s one-third rule. Thus 

$= 0(s) (16) 

(4) From equations (15) and (16) V and $ are known functions of 
s so that 

V = V($) (17) 

Thus V is a known function of $ along the channel wall boundaries in 
the transformed $W-plane. 

(5) If the prescribed velocity distribution along one wall is 
different from that along the other, the channel will, in general, turn 
the flow. This turning angle A0 is given by equation (E5) in 
appendix E. If the turning angle is unsatisfactory a new distribution 
of velocity as a function of s (equations (l4a) and (l4b)) is pre- 
scribed and steps (l) to (5) repeated until the desired value of A0 is 
obtained. Equation (E5) is integrated numerically using Simpson's one- 
third rule, for example ^ .and equation (17) . 

(6) The channel wall boundaries are straight parallel lines of 

constant W equal to 0 and -2£, and extending to in the 

2 , » 
^-direction. Along these boundaries' of constant ¥ a series of equally 
spaced points are located at each of which, the flow direction 6 and the 
physical x_,y-coordinates will he determined hy numerical integration. 

In order to use the tables of a and p presented in this* report; the 
point spacing A<l> must be an even multiple of it/24. Thus the smallest 
point spacing rt/24 is equal to l/l2 of the channel width (n/2) . For 
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a particular prescribed velocity distribution along the channel walls 
the accuracy of the solution increases, and so does the amount of com- 
puting, as the point spacing is reduced. The error for a given point 
spacing- depends on the prescribed velocity distribution and its order 
of magnitude is given by the leading term of the error series of the 
formula used for numerical integration (table VIII, reference 2, for 
example) . For the numerical exan5)le presented in this report the point 
spacing a 4> was jt/l2. From equation (17) 

Alcg^ V (log^ - (log^ V)^ 

A4> AO ' (l®) 

where the subscripts 0 and O+AO refer to adjacent points along the 
channel boundaries, 

(7) The value of 6 at each point channel wall 

A log V 

boundaries is obtained from equation (l3a) in which — is given 

by equation (18) and AL is given by equations (l3b) , (l3c), and table I. 
Note that in equation (l3a) the origin has been moved to 0^ by changing 

A logg V 

from 0 to (O-i^) . Thus the value of for a given value of 

(0-^) varies- with 0^. 

(8) The physical x,y- coordinates at each point on the channel wall 

boundaries are obtained by the numerical integration of equations (5a) 
and (5b) for incompressible flow, or equations (4a) and (4b) for linear- 
ized compressible flow where is given by equation (2e) . con- 

stants of integration in equations (4) and (5) are selected to give known 
values of x and y at upstream or downstream positions where flow 
conditions can be considered uniform. 


NUMERICAL EXAMPLE 

The channel design method of this report has been applied to the 
design of an elbow for the same conditions as exaii 5 )le IV of Part I. 

The design is for an accelerating elbow with no local decelerations of 
the prescribed velocities along the channel -walls and with linearized 
compressible flow. 

Prescribed velocity distribution . - The prescribed velocity dis- 
tribution along the channel -walls is. given by q^ do-wnstream of the 
elbow and by Q as a function of s along the elbow walls. The down- 
stream velocity, q^ is 0.80176. Along the inner wall (with smaller 

radii) of the elbow the arbitrarily prescribed velocity Q as a function 
of arc length s is given by . 
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Q = 0.5 


(s < 0) 




‘ 1 

Q=2+-6"‘^ (0<s< 3.0) 

Q = 1.0 (s > 3.0) 

This velocity distribution is plotted in figure 2. 
From equations (l), (2j), (3c), and (3d) 


J 


(19) 



^2*11 

At* 


Qds 


which together with equation (19) integrates to give 

rt ^2<ld 




2 


(0.5 s) 


(s < 0) 




^ jt ^2^d (b s^ \ ^ ^ r, 

2 "Ap \2 18 ■ loaj (0-S_3.0) 

(s > 3.0) 


'* = f ^ 1-0.75+s) 


J 


( 20 ) 


where from equations (2h), (2i), and (2k) the constant k 2 is equal to 
1.36332 and from equations (2e) to (2k) the constant At* is equal to 
0.73782. From equations (l), (2j), (6c),. (19), and (20) the variation 
in logg V with $ was obtained and is plotted in figure 3. 


The distribution of velocity as a function of arc length is the 
sAme for both channel walls, but, as indicated in figure 3, the dis- 
tribution on the outer wall (larger radii in xy-plane) is shifted in the 

positive $ -direction an amount equal to jt relative to the distri- 

bution on the inner wall. Thus, a velocity difference exists on the two 
walls at equal values of in the interval 0< $5. 3.333n, as shown 
in figure 3. The greater this difference in velocity or the greater the 
range in $ over which velocity differences exist, the greater is the 
elbow turning angle. For the prescribed velocity distribution given in 
figures 2 and 3 the elbow turning angle given by equation (E5) is 
-104.08°. ' 
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Results . - The elbow design resulting from the prescribed velocities 
given in figures 2 and 3 is plotted in figure 4. As indicated in 
table II the contour of this elbow is very nearly the same as that 
obtained by relaxation methods for linearized con^jressible flow with the 
same prescribed conditions (exan^jle IV, Part l) . 

The solution obtained by Green's function (Part II) required one 
experienced computer 3 days whereas the solution by relaxation methods 
(Part I) required about 10 days. The relaxation solutions provide 
additional information, such as the distribution of velocity across the 
channel, but for the most part this additional information is of second- 
ary importance and the design of channels by Green's function is more 
rapid and therefore to be preferred over the design by relaxation 
methods . ' 


SUMMARY OF RESULTS 

Methods of solution by Green's function are developed for the design 
of two-dimensional unbranched channels with prescribed velocities, as a 
function of arc length along the channel walls. The methods apply to 
incompressible and linearized-conpressible, nonviscous irrotational flow. 
One numerical example is presented for an accelerating elbow with 
linearized compressible flow. The elbow shape obtained from the solution 
by Green's function is the same as that obtained from a solution by 
relaxation methods for the same prescribed conditions . The time required 
for the calculations was considerably less for the solution by Green's 
function. 


« 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, September 6, 1951 
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APPENDIX A 
SYMBOLS 

The following symbols are used in this report: 
, . Bernoulli's numbers 


constant, equation (B3) 

Green' s function of the second kind, equations (B2) 
and (10) 

integral (a or p) 

coefficient, equation (2g) 

coefficient, equation (2k) 

length of closed boundary 

distance in xy-plane measured normal to direction of flow 
(expressed as ratio of characteristic length equal to 
channel width downstream at infinity) 

velocity (expressed as ratio of characteristic velocity 
equal to constant channel velocity downstream at 
infinity) ' 

velocity (expressed as ratio of stagnation speed of 
sound) 

velocity used in linearized compressible flow and related 
to q by equation (2J) 

distance from any point in $W-plane to point ($q,Wq) 
at which logarithmic singularity exists 

distance in xy-plane measured along direction of flow 
(expressed as ratio of characteristic length equal to 
channel width downstream at infinity) 

velocity parameter defined by equations (6b) and (6c) for 
incompressible and linearized compressible flow, 
respectively 
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w,Wi,W2 


x,y 


z 

z 


complex functions defined by equations (C3), (Cla), and 
(C2a) , respectively 

Cartesian coordinates in physical plane (ej5)ressed as 
ratios of characteristic length equal to channel width 
downstream at infinity) 

complex coordinate, equation (Clb) 

conjugate, of z 


a 

3 

r 

A 

e 

A0 

p 

p* 

$ 

cp and cp* 
¥ 


¥ and ¥ * 
A¥* 


integral, equation (13d) 
integral, equation (I3e) 
ratio of specific heats 
finite increment 

flow direction in physical xy-plane (measured in t:ounter- 
clockwise direction from positive x-axis) 

channel turning angle, equation (El) 

density (expressed as ratio of stagnation density) 

density used in linearized compressible flow and related 
to p by equation (2f) 

velocity potential used as Cartesian coordinate in trans- 
formed <I>¥-plane and related to cp or cp* 'by equa- 
tion (3a) or (3c) , respectively 

velocity potential for incompressible and linearized 
compressible flow, respectively, equations (3b) and (3d) . 

stream function used as Cartesian coordinate in trans- 
formed $¥-plane and related to ¥ or ¥* ty equa- 
tion (2a) or (2c), respectively 

stream function for incompressible and linearized com- 
pressible flow, respectively, equations (2b) and (2d) 

boundary value of ¥*; for linearized compressible flow, 
along left channel wall when faced in the direction of 
flow, equation (2e) 
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0) 


any harmonic function in $#-plane 


Subscripts: 

a,b quantities related to any two selected values of velocity 

(q^ and q^, respectively) for which densities given 
by equations (2h) and (2i) are equal to densities given 
by equations (2f) and (2g) 


d 

o 

u 




0 


conditions downstream at infinity 

point in $¥-plane at which 0 is determined 

conditions upstream at infinity 

point at on either channel wall boundary 

point at [(‘^*-^o) + either channel wall boundary 

boundary along which W equals 0 


n 

2 


boundary along which ¥ equals 


£ 

2 
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APPENDIX B 


INTEGRAL EQUATION FOR 0($q,¥q) 


If the distribution of the angle in the transformed 

^(If'-plane is harmonic, that is, satisfies equation (8) within and on the 

channel walls (Y equals 0 and ^) , then from Green's theorem and 

the theorem of mean value it can be shown that the value of 0 at .a 
point within (or on) the channel walls is given by (refer- 

ence 3, p. 204, for example) 




2jt 



0^ r ^ 


d$ - 


n 

2 J 


(Bl) 


where the two integrals on the right side of equation (Bl) represent 
the line integral around the channel walls in the counterclockwise 

direction with the signs adjusted so that represents the inner 

o¥ 

normal to the path of integration. 

The function G(iI>,¥) in equation (Bl) is of the form (refer- 
ence 3, p. 204). 


G($,W) = loge i + 03($,¥) 


(B2) 


where r is the distance from any point ($,¥) to the point ($Q,ijfQ) 

and where cd($,¥) is an arbitrary function that is harmonic within and 
on the channel walls. (Thus from equation (B2), G($,¥) is harmonic 

within and on the channel walls except at the point (^o^¥o) where a 
logarithmic singularity exists.) Because the harmonic function 03($,'W) 
is arbitrary, the function G(<I>,¥) can be selected so that along the 

tr 

channel wall boimdaries (i^r equals .0 and is a constant c 

^ o¥ 

given by the following equation (obtained from notes presented by 
Tamarkin and Feller in the 1941 Summer Session for Advanced Instruc- 
tion and Research in Mechanics at Brown Univ.): 


2 « 

I 


(B3) 


where I is the length of the path along which the line integral is 
taken. For the path under consideration I is infinite and therefore 
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can "be selected so that is zero along the channel walls. 

A function with this property is called a Green' s ftinction of the sec- 
ond kind. Equation (Bl) becomes 


®(<5o^^o) 


or, combined with equation (6a) 







0 


d4» 




Along the channel walls 



(“ - (“ 


d$ 


( 9 ) 


a log^ V 




— is known from the prescribed veloc- 


ity distribution so that, after the proper Green's f-unction G has 
been determined (appendix C), equation (9) determines the value of 0 
at any point ($Q,iirQ) . The value of 0($Q,fQ) given by equation (9) 

can be adjusted by an arbitrary constant of integration to give a spec- 
cified value of 0 at one point in the flow field. 
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APPENDIX C 

GEEEN'S FUNCTION OF SECOND KIND 

From appendix B Green's function of the second kind G satisfies 
the condition 


^ = 0 


along the channel walls, which are straight and parallel boundaries 
(iy equals 0 and extending to ±oo in the $ -direction, and 
satisfies the equation 




everywhere in the channel except at the point G has 

a logarithmic pole. For these conditions the Green’s function G can 
he obtained by analogy from the velocity potential for incompressible 
flow into a point sink at ’ between straight parallel bound- 
aries at equal to 0 and The logarithmic pole for G at 

(®o^^o) corresponds to the point sink and the condition — a 0 at 

oV 

the boundaries corresponds to zero velocity, that is, no flow normal 
to the boundaries. 


The velocity potential for fluid flow with the boundary conditions 
just described is obtained from two infinite series of point sinks with 
the sinks of each series spaced n distance apart in the ijr-direction 
and the two series arranged by the method of images in such a manner 

that no flow crosses the boundaries, . that is ^ = 0. This arrangement 
of point sinks is shown in figure 5. 

The complex function wj^ for the first infinite series of point 
sinks is given by (reference 4^ p. 112, for example) 

w^ =s -logg sinh (z-Zq) (Cla) 

where 


z == $ + iW 


(Clb) 
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The complex function ¥3 for the second infinite series of point sinks 
(mirror image of the first series in order to prevent flow across the 
boundaries W equals 0 and is given by 

Wg = -logg sinh (z-Zq) (C2a) 

where 

z = ^> - iV (C2b) 

The complex function w for the combined flow becomes from equa- 
tions (Cla) to (C2b) , , 

w = w^ + Wg = -logg sinh [($-^> 0 ) + i " 

logg sinh [(0>-$g) + i (¥+¥g)] (C3) 

The Green's function of the second kind G corresponds to the 
velocity potential for the incompressible flow and is therefore given by 
the real part of equation (C3) 

G = - loge [cosh^ ($-<I>g) - cos^ (^-■ 3 ^ 0 )] [cosh^ (^-^o) - (Mg)] 

(C4) 

But along the channel walls ¥ is equal to 0 or ^ so that 

”V 

cos^ (¥+¥g) = cos2 (¥-¥g) 

and equation (C4) becomes 

Go or £ = -loggfcosh^ ($-$o) - cos^ (¥-¥g)] (lO) 

2 - 

Equation (lO) gives the Green's function of the second kind along the 
channel walls (straight parallel lines of constant ¥ equal to 0 and 

— and extending to ±<» in the ^-direction) . 
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APPENDIX D 


EVALUATION OF a AND P 

Several techniques^ depending on the magnitude of the upper limit 
, were used to evaluate the integrals a and p given by- 
equations (I3d) and (l3e) . Each integral is treated separately in this 
appendix and the values of for the upper limit |(^^-^b)l are 

considered positive. For negative values of ($-$ 0 ) magnitudes 

of I (that iS; of a or p) are equal for corresponding values of 
|($-$q)| but opposite in sign. As a result the values of AI have the 
same sign. 


Integral a 

Small and medium values of ('^*-'^ 0 ) • - For small and medium values 
of the -upper limit of integration ($- $q) in equation (l3d), that is, 
for 0 < < 60 jt/24, the integral a is evaluated by Simpson's 

one-third rule using increments of equal to rt/48. 

Large values of . - For large values of (®-%) that is 

for ($- <^) > 60 Jt/24, the integrand in equation (l3d) becomes 

logg cosh - logg 2 (Dl) 


so that equation (l3d) becomes 
\ort/24 


a s 


logg cosh d(5>-'J’g) + 


^60jt/24 


- log^ 2] d($-$^) 


» 25.809782 + 


[ 




- 0.693147 - 25.398552 




« 0.411230 - 0.693147 ($- $ ) + i ($-<!> f 


(D2) 
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Equation (D2) gives values of a for values of equal to or 

greater than 60n/24. Values of the integral a are tabulated in table I 
for a range of |($-$q)| between 0 and lOOit/24 in increments of jt/24. 

For negative values of the sign of a is negative. 

Integral p 

SmaOJ. values of • ~ For (O-^q) equal to zero the integrand 

of equation (13e) becomes minus infinity so that Simpson's one-third 
rule cannot be used to evaluate p in this region of , as was 

done for a. However^ equation (I3e) integrates by parts to give 

logg sinh ($-$q) d($-$g) 

= logg sinh - 




(0-$g) ctnh d($-^) 


(D3) 


where the integrand ctnh on the right side of equa- 

tion (D3) can be expanded in the following series form: 


22b, 

($-<I)o) ctnh = 1 + O' 


2^B3 S^Bg 

H 61 


2®B., (<!>-$ )® 210Bq 212b.^ ($-$)12. 

Z_ + L_ 2L- iL_ + (m) 

81 ^ 10! 121 + . . . 

where B^^, B3, and so forth, are called Bernoulli's numbers (refer- 
ence 5, p. 90, for example). From equations (D3) and (D4) 
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p = (<I>-^o) logg sinh - i^-%) - g + 


225 


2($.$^)7 
6615 ■*■ 


1382 

42,525 1,029,105 ■*■ 8,300,667,375 " • * * (®5) 


Equation (D5) was used to obtain p .as a function of . for 

0 < < 8jt/24. 


tfedium values of (3 >-$q) . - For medi um values of the upper limit 

of integration in equation (l3e), that is, for 8jt/24 < . ($-0^) 

< 60rt/24, the integral p is evaluated by Simpson' s one- third rule as 
was done for a. 


Large values of • - For large values of that is, 

for > 60rt/24, the integrand in equation (l3e) becomes 


logg sinh » ($-$o) - logg 2 

so that equation (l3e) becomes 
'60rt/24 


(D6) 


P s 


'0 


log^ sinh (0-$^) d($-^) + 


60jt/24 


[($-$o) - loSe 2] d(<^-$o) 


*» 24.576082 


+ |_ 2 0.693147 ($-$o) - 25.398552J 


0.822470 - 0.693147 | (D7) 

Eqiiation (D7) gives values of ^ for values of (<1’-$^) equal to or 

greater than 60jt/24. Values of the integral p are tabulated in table I 
for a range of |('5 -$q)| between 0 and lOOn/24 in increments of n/24. 

For negative values of ($-$q) the sign of ^ changes. 
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APPENDIX E 


CHAKHEL TURNING ANGLE 


If the prescribed velocity distribution along one channel wall 
differs from the distribution along the other wall, then in general the 
channel deflects an amount A0, which is the difference in flow direc- 
tion far downstream and far upstream of the region in which the pre- 
scribed velocity distribution varies. Thus, 


A0 = 0a - ®u 

For large values of' ICt-fo)! such as occur far upstream and far 

downstream of the region in which the prescribed velocity varies along 
the channel walls 

cosh^ » cos2 


so that from equation (lO) 


^ ^ 

2 


p(<K-$o)| - logg 2j 


(E2) 


Far upstream that 


|(4.-$o) I = (*-»„) 

and because V is harmonic 



so that equation (E2) substituted in equation (9) gives 




( S loge V' 


E 

2 


/ h loge V \ 

\ 5 $^ Jo 


d$ 


(E3) 
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Likewise, far dovmstream so that 




and equation (E2) substituted in equation (9) gives 


Od 




From equations (El), (E3), 


and (E4) 



Equation (E5) determines the channel turning angle A0. 


(E4) 


(E5) 
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TABLE I - TABULATED VALUES OF THE INTEGRALS a AND P FOR A RANGE OF |(<I'-<Ib)l 
[Computational methods given in appendix D.J 



CL 

0 

0 

l{«/24) 

0.000373 

2{it/24) 

.002970 

3(n/24) 

,009942 

4(jt/24) 

.023301 

5(7t/24) 

.044875 

6(k/24) 

.076277 

7(«/24) 

.118397 

8(n/24) 

. 173909 

9(jt/24) 

.242283 

10(n/24) 

,324812 

ll(«/24) 

.422136 

12(rt/24) 

.534763 

13(n/24) 

. 663098 

14{jt/24) 

.807460 

15(ji/24) 

.968096 

16(«/24) 

1, 145201 

17(n/24) 

1.338926 

18(n/24) 

1.549389 

19(it/24) 

1.776679 

20(«/24) 

2.020867 

21(jt/24) 

2.282008 

22(«/24) 

2.560143 

23(n/24) 

2.855304 

24(Tt/24) 

3.167516 

25(rt/24) 

3.496799 

26(n/24) 

3.843168 

27(jt/24) 

4.206633 

28(n/24) 

'4.587203 

29{n/24) 

4.984885 

30(rt/24) 

5.399685 

31(rt/24) 

5.831607 

32(i(/24) 

6.280652 

33(jt/24) 

6 . 746825 

34(jt/24) 

7.230126 

35(ir/24) 

7 , 730557 

36(«/24) 

B. 248119 

37(«/24) 

8.782813 

38(it/24) 

9.334640 

39(it/24) 

9.903600 

40(ji/24) 

10.489693 

41(n/24) 

11.092920 

42(it/24) 

11.713281 

43(n/24) 

12.350777 

44(jt/24) 

13.005406 

45(n/24) 

13.677170 

46{jt/24) 

14.366068 

47(n/24) 

15.072101 

48(*/24) 

15.795268 

49(fl/24) 

16.535571 




(^«/ 24 ; 


)|(44>=2»/Z4) 


.013359 


.042620 


.^46569 


.451922 


.209951 


.337741 


1.360662 


fl(l) 


-.954346 


-.936430 


- . 755399 


.792488 


1.935683 


2.973783 


4.598205 


5.047181 


6.496961 


8.669936 


11.771715 


13.132373 


AI = A3 


A^P 

(Al»=n/24) 

AgP 

(A<N2n/24) 

-0.396937 

-0.611660 

-.214723 

- . 144746 

-.242791 

-.098045 

-.062035 

-.094079 

- . 032044 

-.005816 

.012092 

.017908 

.039896 

. 100542 

.060646 

.080497 

.180181 

.099684 

.118377 

.255075 

. 136698 

. 154739 

.327305 

.172566 

. 190232 

.398006 

.207774 

.225221 

,467817 

.242596 

.259915 

.537106 

.277191 

.294435 

.606089 

.311654 

.328853 

.674890 

. 346037 

. 363210 

. 743585 

.380375 

.397531 

.812215 

, ,414684 

.431832 

.880808 

.448976 

.466120 

.949380 

.483260 

.500400 

1.017938 

.517538 

.534676 

1.086488 

.551812 

■ .568949 

1.155034 

■ .586085 

■ .603220 

1.223576 

■ .620356 

.637492 

1.292118 

■ .654626 

■ .671762 

1.360658 

' . 688896 

“ .706032 

1.429198 

.723166 

■ . 740302 

1.497738 

“ .757436 


( 1 ) 


For negative values of (4-®o) "tbe signs of a and 3 change, hut the signs 
of Aa and A3 remain unchanged. 
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TABLE I - TABULATED VALUES OF THE IBTEGRALS a AHD S FOR A RAHGE OF |($-*g)| - Concluded 
[Computational methods given In appendix D.3 



a 

50(«/24) 

17.293007 

51(jr/24) 

18.067579 

52(jt/24) 

18.859285 

53(it/24) 

19.668125 

54( It/24) 

20.494101 

55(jt/24) 

21.337211 

56{it/24) 

22.197456 

57(jt/24) 

23.074835 

58(jt/24) 

23.969349 

59 ( It/24) 

24.880998 

60(it/24) 

25.809782 

6l(it/24) 

26,755694 

62(it/24) 

27.718750 

63(it/24) 

28.698940 

64(jt/24) 

29.696257 

65(k/24) 

30.710717 

66(it/24) 

31.742311 

67(it/24) 

32.791033 

68(it/24) 

33.856896 

69(it/24) 

34.939895 

70(it/24) 

36.040029 

71(it/24) 

37.157289 

72(it/24) 

38.291692 

73(jt/24) 

39.443230 

74(it/24) 

40.611893 

75(it/24) 

41.797701 

76(it/24) 

43.000643 

77(it/24) 

44.220710 

78{it/24) 

45.457922 

79(it/24) 

46.712269 

80(it/24) 

47.983750 

81(it/24) 

49.272356 

82(it/24) 

50.578107 

83( It/24) 

51.900992 ■ 

84{it/24) 

53.241002 - 

85(it/24) 

54.598158 - 

86(it/24) 

' 55.972447 ■ 

87(jt/24) 

57.363861 - 

88{it/24) 

58.772421 - 

89(it/24) 

60.198115 - 

90(it/24) 

61.640933 - 

91(it/24) 

63.100897 ' 

92(it/24) 

64.577995 - 

93(rt/24) 

66.072228 ■ 

94(it/24) 

67.583585 - 

95(n/24) 

69.112088 - 

96(jt/24) 

70.657725 ~ 

97(jt/24) 

72.220486 _ 

98(it/24) 

73.800393 - 

99 ( It/24) 

75.397435 - 

100(it/24)^2^ 

77.011599 ■ 


(A<ft=rt/24)| (A4fc2«/24) 


.774572 


.791706 


,808840 


.825976 


,843110 


,860245 


,963056 


.980190 


1.014460 


1.065863 


1.100154 


1,117260 


1.134403 


1.151538 


1.168663 


1.185808 


1.220067 


1.254347 


1. '442818 


1.459964 


1.477098 


1.528503 


1.545637 


1.579907 


1.597042 


1.614164 


1.771893 


1.840433 


1.977507 


2.046044 


2.114585 


2.183123 


2.251663 


2.388740 


2.525818 


3.005590 


3.211206 


p( i) 

16.059309 

16.833880 

17.625585 

18.434426 

19.260401 

20.103511 

20.963755 

21.84113.6 


22.735649 

23.647298 

24.576082 

25.521994 

26.485050 

27.465240 

28.462557 

29.477017 

30.508611 

31.557333 

32.623196 

33.706195 

34.806329 

35.923589 

37.057992 

38.209530 

39.378193 

40.564001 

41.766943 

42.987010 

44.224222 

45.478569 

46.750050 

48.038656 

49 . 344407 . 

50.667292 - 

52.007302 . 

53.364458 - 

54.738747 - 

56 . 130161 - 


57.538720 - 

58.964415 - 

60.407233 - 

61.867197 - 

63.344295 - 

64.838528 - 

66.349885 - 

67.878388 - 

69.424025 - 

70.986786 

72.566693 - 

74.163735 - 

75.777899 ~ 


AI = A6 


AiP 

(A^jt/24) 


.791705 


.877380 


.980190 


.997517 


1.082999 


1.134403 


1.151538 


1.166665 


1.185808 


1.202942 


1.220067 


1.237212 


1.254347 


1.271481 


1.408560 


1.494233 


1.579907 


1.597042 


1.614164 


( 1 ) 


( 2 ) 


For negative values of the signs of 

of Aa and A8 remain unchanged. 


(A^2«/24) 


1.840433 


1,908968 


1.977507 


2.046044 


2.114585 


2.183123 


2.251663 


2-320201 


2.457279 


2.662895 


2.937052 


3.074130 


3.142668 


For values of 
for p. 


and p change j hut the signs 
>100(it/24) use equation (D2) for a and equation (D7) 
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TABLE II - C(»IPABISON OF ELBOW IfBSKKS OHTAIHED FROM SOLUTIONS BY RELAXATICHI METHODS AND BY QIEEN'S FUNCTION 


[Linearized compressiLle flow; prescribed velocity distribution given in figs. 2 and 3] 






y = 0 

(inner wall) 


¥ = it/2 
(Outer wall) 

Q 

q 

Tela 

olutlon by 
xation methods 
(Part I) 

Solution by 
Green's function 
(Part II) 

Q 

q 

Solution by 
relaxation methods 
(Part I) 

Solution by 
Green's function 
(Part II) 

X 

y 

e 

(deg) 

X 

y 

e 

(deg) 

X 

y 

e 

(deg) 

X 

y 

e 

(deg) 

-22(rt/24) 

0.5000 

0.4009 

-2.466 

-0.769 

0 

-2.466 

-0.769 

0 

0.5000 

0.4009 

-2.466 

.770 

0 

-2.466 

.770 

0 

-20(n/24) 

.5000 

.4009 

-2.241 

-.769 

.01 

-2.241 

-.769 

.01 

.5000 

.4009 

-2.241 

.770 

-.01 

-2.241 

.770 

-.01 

-18( It/24) 

.5000 

.4009 

-2.016 

-.769 

.01 

-2.016 

-.769 

.01 

.5000 

.4009 

-2.016 

.770 

-.01 

-2.016 

.770 

-.01 

-16U/24) 

.5000 

.4009 

-1.791 

-.769 

.02 

-1.791 

-.769 

.02 

.5000 

.4009 

-1.791 

.770 

-.02 

-1.791 

.770 

-.02 

-14(n/24) 

.5000 

.4009 

-1.566 

-.768 

.03 

-1.566 

-.768 

.03 

.5000 

.4009 

-1.566 

.770 

-.03 

-1.566 

.770 

-.03 

-12(«/24) 

.5000 

.4009 

-1.341 

-.768 

.05 

-1.341 

-.768 

.05 

.5000 

.4009 

-1.^41 

.769 

-.05 

-1.341 

.770 

-.05 

-10(«/24) 

.5000 

.4009 

-1.116 

-.768 

.08 

-1.116 

-.768 

.08 

.5000 

.4009 

-1.116 

.769 

-.08 

-1.116 

.769 

-.08 

-8(n/24) 

.5000 

.4009 

-.891 

-.768 

.14 

-.891 

-.768 

.14 

.5000 

.4009 

-.891 

.769 

-.13 

-.891 

.769 

-.13 

-6( It/24) 

.5000 

.4009 

-.666 

-.767 

.24 

-.666 

-.767 

.24 

.5000 

.4009 

-.666 

.768 

-.21 

-.666 

.768 

-.22 

-4(«/24) 

.5000 

.4009 

-.441 

-.766 

.40 

-.441 

-.766 

.41 

.5000 

.4009 

-.441 

.767 

-.35 

-.441 

.767 

-.36 

-2(«/24) 

.5000 

.4009 

-.216 

-.763 

.70‘ 

-.216 

-.763 

.74 

.5000 

.4009 

-.216 

.765 

-.56 

-.216 

.765 

-.59 

0 

.5000 

.4009 

.008 

-.760 

1.31 

.010 

-.759 

1.52 

.5000 

.4009 

.009 

.763 

-.92 

.009 

.762 

-.94 

2(it/24) 

.5079 

.4072 

.233 

• -.752 

. 2.82 

.233 

-.751 

2.76 

.5000 

.4009 

.234 

.758 

-1.45 

.234 

.758 

-1^48 

4(n/24) 

.5293 

.4243 

.450 

-.739 

3.88 

.449 

-.738 

3.76 

.5000 

.4009 

.459 

.751 

-2.22 

.459 

.750 

-2.24 

6{ It/24) 

.5599 

.4489 

.656 

-.724 

4.00 

.656 

* -.724 

3.89 

.5000 

.4009 

.684 

.740 

-3.28 

.684 

.740 

-3.29 

8(«/24) 

.5962 

.4780 

.851 

-.712 

3.17 

.850 

-.713 

3.07 

.5000 

.4009 

.908 

.725 

-4.65 

,908 

.724 

-4.68 

10( It/24) 

.6354 

.5094 

1.033 

-.704 

1.44 

1.033 

-.705 

1.38 

.5000 

.4009 

1.132 

.703 

-6.41 

1.132 

.702 

-6.43 

12(it/24) 

.6754 

.5415 

1.205 

-.703 

-.98 

1,205 

-.704 

-1.04 

,5000 

.4009 

1.355 

.674 

-8.52 

1.355 

.673 

-8.56 

14(it/24) 

.7149 

.5732 

1.366 

-.710 

-4.00 

1.367 

-.711 

-4.04 

.5000 

.4009 

1.577 

.636 

-11.04 

l'.577 

.635 

-11.06 

16(n/24) 

.7531 

.6038 

1.519 

-.725 

-7.48 

1.519 

-.727 

-7.51 

.5000 

.4009 

1.797 

.587 

-13.91. 

1.797 

.586 

-13.94 

18(it/24) 

.7894 

.6329 

1.663 

-.749 

-11.35 

1.663 

-.750 

-11.37 

.5000 

.4009 

2.013 

.527 

-17.13 

2.013 

.526 

-17.16 

20(it/24) 

.8235 

.6602 

1.798 

-.781 

-15.52 

1.799 

-.783 

-15.54 

.5000 

.4009 

2.226 

.455 

-20.67 

2.226 

.453 

-20.69 

22(it/24) 

.8550 

.6855 

1.926 

-.822 

-19.96 

1.926 

-.823 

-19.98 

.5000 

.4009 

2.434 
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Figure 3. - Variation in prescribed values of log^ V with ^ along channel walls of numerical example 









• Figure 5. - Two infinite series of point sinks required in the 
development of Green's function of the second kind G. 


NACA -Langley - 1-4-52 - 1000 



